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We consider multi-flavor QED on a finite lattice at a finite chemical potential and show that 
the partition function only depends on the variables, — ^j, for i = 2, ■ ■ ■ , where qi, i — 1 ■ • ■ 
are integer valued charges of the various flavors and [ii , i = 1 ■ • • are the dimensionless chemical 
potentials of the various flavors. 
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Consider multi-flavor QED regularized on the lattice. We will assume that the continuum theory is defined by first 
defining the theory on a finite periodic lattice and then taking the limit where the extent in all directions goes off to 
infinity. We will not be concerned in this paper as to whether such a theory has a well defined continuum limit. For 
definiteness, we could think about two dimensional multi- flavor QED which has a well defined continuum limit. 

Let J7i/(n) be the link variable connecting the sites n and n + v on a d dimensional periodic lattice, L-y x • • • Ld- 

Consider the class of gauge fields given by C/„(n)e L » with < h v < 1 for v = 1, • • • d. Gauge fields with different 
choices of the dimensionless variables, h v , within this class are not gauge equivalent but have the same gauge action. 
The fermion determinant, on the other hand, depends on the variables, h v , which we will refer to as the toron 
variables. Consider a fermion with integer charge qi and a dimensionless chemical potential /i^. The chemical potential 

is introduced pQ by multiplying the parallel transporter in the forward v = d direction by e L <i and in the backward 

v = d direction by e L d . The factor of 27r is introduced for convenience and we keep the dimensionless chemical 
potential, fii, fixed as we take Ld — > oo. 

The fermionic determinant is a function of U v {n), v = 1, • • • ,d; h u for v = 1, - ■ ■ ,d— 1; and z, = hd — i^r- Since the 
gauge action does not depend on the toron variables, we can integrate the fermion determinant over these variables 
using the uniform measure. Consider, for simplicity, a lattice fermion operator of the naive, Wilson or staggered type. 
We will show later that our arguments will also apply to the overlap Dirac operator under certain mild assumptions. 
The lattice fermion operator is a finite matrix on a finite lattice. Focussing on the dependence on hd alone, we see 

that the fermion determinant is a finite polynomial in e L * and e L d . The fermion determinant will be an 
analytic function of Zj in the complex plane. Since the fermion determinant is gauge invariant, it will be periodic 
under Zj — > Zi + 1. A contour integral in the complex plane results in the integral over hd in the range [0, 1] to be 
independent of /i;. If we consider a theory with many flavors, that all have the same value for ^f, then again the 
integral over hd will yield a result that is independent of all the fii . In other words, a multi-flavor theory of QED at a 
finite chemical potential can only depend on the variables, — , for i = 2, • • • . This result will also hold for an 

anomaly free chiral QED as long as the lattice formulation is gauge invariant, at least in the continuum limit. This 
is the main observation in this paper. 

Massless Schwinger model in the presence of a chemical potential was first studied in [2] . The problem was treated 
in the Hamiltonian formalism. In order to deal with a finite problem, a uniform charge background was introduced 
in a finite region of space. This causes an explicit breaking of translational invariance. The ground state is a classical 
Wigner crystal which is not destroyed by quantum fluctuations. An explicit chemical potential term was introduced 
in the Hamiltonian formalism in [3]. Their formalism also had to break translational invariance and they conclude 
that there is a inhomogeneous chiral condensate in the Schwinger model at finite density. A path integral formulation 
of the problem, again with the introduction of a chemical potential that breaks translational invariance results in 
an inhomogenous chiral condensate [4]. References to a possible inhomogenous chiral condensate in the Schwinger 
model are still being discussed in the literature [SH7]. The problem is discussed in [5] where the author argues that 
the inhomogeneous chiral condensate in the Schwinger model at finite density is an artifact of the explicit breaking 
of translational invariance in the formalism. The generalized Thirring model was analyzed in [3] and it has been 
explicitly shown that the chiral condensate does not depend on the chemical potential. Our observation in this paper 
clearly shows that physics does not depend on the chemical potential in the Schwinger model. 
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FIG. 1: The left top and bottom panels show the behavior of the determinant of the overlap Dirac operator in the zero and 
unit topological sectors. The right top and bottom panels show the behavior of the determinant of the Wilson Dirac operator 
in the zero and unit topological sectors. The results are on a 7 x 7 lattice and the dimensionless chemical potential is set to 
and the x axis shows the value of /i2- 

We now proceed to discuss the case of overlap fermions in some detail when the chemical potential is introduced 
as in [10 . Since the definition involves the sign function of a complex matrix, it is not apriori clear if the arguments 
presented above will apply. This problem is addressed using a specific numerical example, namely, the Schwingcr 
model with a finite chemical potential, /i. Gauge fields are generated at a fixed coupling using the gauge action 
described in [11] with zero and unit topological charge. The determinant of both the Wilson-Dirac operator and 
the overlap-Dirac operator was computed in a fixed gauge field background that has zero topological charge. The 
determinant of both the Wilson-Dirac operator and the overlap-Dirac operator (we exclude the zero mode in this 
case) was also computed in a fixed gauge field background that has unit topological charge. Averages of these two 
quantities over all gauge fields for the overlap-Dirac operator enter the computation of the chiral condensate in the 
masslcss Schwinger model. 

We first consider two sample gauge fields, one with zero topological charge and one with unit topological charge as 
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FIG. 2: The left top and bottom panels show the behavior of the determinant of the overlap Dirac operator in the zero and 
unit topological sectors. The right top and bottom panels show the behavior of the determinant of the Wilson Dirac operator 
in the zero and unit topological sectors. The results are on a 7 x 7 lattice and the dimensionless chemical potential is set to 
0.5 and the x axis shows the value of hi. 

set by the gauge fields, on a 7 x 7 lattice. The same gauge field is used for several different values of the chemical 
potential. Plots of the real part of the determinants [17] for three sample values of the chemical potential, namely, 
fj, = 0,0.5, 1.0 arc shown in Figs ] 1|2|3] respectively. In all three cases, the determinant of the Wilson Dirac operator 
with a specific choice for the mass term used in the definition of the overlap Dirac operator kernel behaves smoothly 
as a function of hi. The determinants do depend on hi and the variation increases as one increases the chemical 
potential. Furthermore, the function even changes sign for /i = 1.0. In spite of this, we explicitly verified that the 
integral over hi (which remains complex for a fixed gauge field background) is independent of \x as expected by the 
analytical argument presented in the beginning of the paper. 

The sign function of the Wilson Dirac operator need not be a smooth function of hi. This is due to the fact it 
depends on the sign of the real part of the eigenvalue of the Wilson Dirac operator and this can change as a function of 
hi . The trace of the sign function of the Wilson Dirac operator is defined as twice the topological charge of the gauge 



4 




'L.mtLy 

0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



FIG. 3: The left top and bottom panels show the behavior of the determinant of the overlap Dirac operator in the zero and 
unit topological sectors with zero chemical potential. The middle top and bottom panels show the behavior of the determinant 
of the Wilson Dirac operator in the zero and unit topological sectors. The right top and bottom panels shows the topological 
charge measured by the overlap Dirac operator in the presence of a chemical potential for the presumed zero and unit topological 
charge configuration. The results are on a 7 x 7 lattice and the dimensionless chemical potential is set to 1.0 and the x axis 
shows the value of /12 ■ 

field. Whereas, the results are consistent with the topological charge of the constructed gauge field background for 
small values of the chemical potential, it does not remain consistent for large values of the chemical potential as seen 
in the right panel of Fig. [3] The top panel should have been consistent with zero and the bottom panel should have 
been consistent with unity. Since this is not the case for \i = 1, the determinant of the overlap Dirac operator in the 
two different topological sectors is not a smooth function of hi as seen in the left panels of Fig. [3j As a consequence 
the independence of the integral of the determinant of the overlap Dirac operator on the chemical potential breaks 
down beyond a certain value of the chemical potential. This is due to the discretization of the dimensional chemical 
potential, /i, over a finite number of slices, L^, taken to be 7 in Fig. [3j If we increase the number of slices, Ld, 
keeping the physical gauge coupling fixed, we found that the value of the chemical potential where the independence 
breakdown increases. For example, we could only go up to li — 0.7 on 7 x 7 lattice at a given coupling but we could 
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FIG. 4: The left top and bottom panels show the behavior of the determinant of the overlap Dirac operator in the zero and 
unit topological sectors with zero chemical potential. The right top and bottom panels show the behavior of the determinant 
of the Wilson Dirac operator in the zero and unit topological sectors. The results are on a 9 x 9 lattice and the dimensionless 
chemical potential is set to 1.0 and the x axis shows the value of /i2- 

go up to fi = 1.0 on a 9 x 9 lattice at the same physical coupling. This is illustrated for the case of fi = 1.0 on a 
9x9 lattice where the behavior of the determinants of the overlap Dirac operator are smooth as a function of /12. 
The problem with the overlap Dirac operator in the presence of a chemical potential was anticipated in |10j but is 
expected not to affect the continuum limit. 

We have presented an analytical argument in this paper that there is no dependence on the chemical potential in 
the path integral formalism of QED. This result is valid within the lattice formulation and in the continuum limit. 
Whereas, there is some justification to the argument presented in [5] toward the problems with breaking translation 
invariance in the Hamiltonian formalism, we have shown that the main reason for recovering the correct behavior 
in the presence of a chemical potential is the integration over the toron variable as emphasized in [3] and further 
emphasized in |12j . One does not see the toron variable, hd, in the hamiltonian formalism since one starts in the 
Coloumb gauge. The realization of the toron variable in the Hamiltonian formalism would be to integrate over all 
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boundary conditions for fermions in the Euclidean time (temperature) direction |13j . 

We have not performed an analytic continuation in the chemical potential in this paper. One could reproduce the 
central argument in this paper by working with an imaginary chemical potential |18j . Since there is no periodicity 
when the chemical potential is real, we refrained from using imaginary chemical potential - our results are valid for 
all values of the real chemical potential. A non-trivial dependence on the chemical potential will be seen in the two 
flavor Schwinger model and finite density phase transitions are expected in four dimensional two flavor QED. 

The continuous toron variable becomes a discrete Zjq variable in a SU(N) gauge theory and our argument of 
independence on the chemical potential will not go through [19] In the limit of N — > oo, we have a continuous toron 
variable. Therefore, we expect the toron variable to play a part in the analysis of the 't Hooft model in the presence 
of a chemical potential as discussed in [16] . The Gross-Neveu model [6] is different in this aspect since one does 
not integrate over all possible fermionic boundary conditions in the Euclidean time direction. This will be the case 
even if we introduce a bosonic variable to convert the four-fermi coupling into a fermion bilinear since we will have a 
Gaussian term for the bosonic field. 
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If we analytically continue to imaginary chemical potentials, the issue of Roberge- Weiss transitions [15] need to be adressed. 



